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Chameleon mechanism appearing in massive tensor-scalar theory of gravity can effectively reduce 
locally the non-minimal coupling between the scalar field and matter. This mechanism is invoked 
to reconcile large-scale departures from general relativity (GR), supposedly accounting for cosmic 
acceleration, to small scales stringent constraints on GR. In this paper, we carefully investigate 
this framework on cosmological and solar system scales to derive combined constraints on model 
parameters, notably by performing a non-ambiguous derivation of observables like luminosity dis- 
tance and local post-newtonian parameters. Likelihood analysis of type la SuperNovae (SN la) data 
and of admissible domain for the PPN parameters clearly demonstrates that chameleon mechanism 
cannot occur in the same region of parameters space than the one necessary to account for cosmic 
acceleration with the assumed Ratra-Peebles potential and exponential coupling function. 
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I. INTRODUCTION 

Today, gravitation is facing a major problem: on one 
hand. General Relativity (GR) have passed all solar 
system experiments ; on the other hand, GR and the 
Standard Model of particles are not sufficient to explain 
galactic and cosmological observations, for instance the 
accelerated expansion of the universe. Several options 
have been advanced to solve this apparent inconsistency. 
The most widespread possibility consists of extending 
the matter-energy content of the universe with new 
ingredients like dark matter and dark energy. Another 
possibility consists of modifying GR on large scales 
without the introduction of a new type of matter. The 
construction of a new theory of gravity is particularly 
hard because it has to satisfy very stringent constraints 
on solar system scales (see Will [I] for a review) while 
showing deviations on large scales. 

Moreover from a theoretical point of view, the different 
attempts to quantize gravity or to unify it with other 
fundamental interactions predict deviations from GR. 
In most of the theories proposed to account for cosmic 
acceleration, to quantize gravity or to unify it with 
the standard model, scalar fields naturally appear in 
addition to the metric tensor. For examples, scalar fields 
appear in Kaluza-Klein theories or in string theories. 
They also appear in non-minimal extension of GR like 
in f{R) theories [5] or in credible cosmological scenarios 



(infiation or quintessence). 

Tensor-scalar theories are therefore very important 
and they are already widely studied in the literature. 
The first tensor-scalar theory was introduced by Jordan, 
Brans and Dicke [3', '4^ in order to recover the Mach princi- 
ple. A very detailed study of tensor-scalar theory can be 
found in Damour and Esposito-Farese [S]. Tensor-scalar 
theories incorporate both the direct coupling of the scalar 
field to matter (through the coupling function A{(j)) and 
a scalar self-interaction term (represented by a potential 
F((/))). Massless tensor-scalar theories (characterized by 
a vanishing potential V{(j)) = 0) are strongly constrained 
at the linear level of the coupling function by solar sys- 
tem experiments [T| and by binary pulsars ^Gj , yet rather 
poorly at the non-linear level. In particular, the observa- 
tion of the Shapiro delay with Cassini spacecraft during 
a solar conjunction [7 gives currently the best constraint 
on the post-newtonian parameter 7 

7- 1 = 2.1 ±2.3 X 10"^- (1) 

This constraint implies that the coupling function at the 
linear level dlog{A)/d4> = k{(j)) < 1 Q This low limit 
puts a severe constraint on some theoretical development 
that requires a coupling of the order of the unity (for 
example string theory or f{R) gravity produce coupling 
of the order of the unity) unless some mechanism is 
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reducing this coupling constant. 

In 2004, Khoury and Weltman |5] showed that 
certain massive tensor-scalar theory can satisfy solar 
system experiments even for high coupling constants. 
This happens because the amplitude of the mass of 
the scalar field depends on the local matter density. 
In particular, in a region of high density (for example 
in the Sun) the scalar field acquires a high mass while 
on cosmological scales (where the local density is low) 
the mass of the scalar field is very low. The fact that 
the scalar field has a mass mimicking the local density 
justifies the name given to these scalar field called 
chameleon fields. 

Considering a spherical body, Khoury and Welt- 
man [HI [5] showed that under certain conditions, the 
scalar field is nearly frozen inside the body at the 
exception of a thin shell around the surface's body. This 
screening effect appears only if the so-called thin-shell 
parameter (noted e) is small enough. This effect has the 
interesting property to reduce the effective scalar charge 
of the body that is equal to the coupling constant k 
in usual tensor-scalar theory in low gravitational field 
(non perturbative effects can appear in strong field as 
shown in Damour and Esposito-Farese [T^). Since the 
effective scalar charge entering in the expression of the 7 
post-newtonian parameters is reduced, this mechanism 
allows the theory to pass the Cassini constraint (and 
more generally solar system experiments) even for high 
coupling constants k that were previously excluded by 
PPN constraints. 

The main constraint on chameleon model on solar 
system scales is that the mass of the scalar field has to 
be sufficiently large to satisfy solar system constraints 
(equivalence principle, fifth force and PPN constraints). 
In Khoury and Weltman '8\, it is explicitly argued that 
with this mechanism it is possible to explain cosmological 
effects without producing observable deviations in solar 
system because the environment is dense enough and 
the contribution from the scalar field is exponentially 
suppressed. 

The goal of this paper is to verify this assumption and 
to derive acceptable domain for the parameters charac- 
terizing the theory (mainly the coupling constant k and 
the parameters characterizing the potential V{(f>)). The 
chameleon mechanism has already been widely studied in 
the literature [8j ^ ITT| - fT5] . Nevertheless some issues need 
still to be addressed. First of all, none of these articles 
derive constraints from Supernovae data taking correctly 
into account non minimal coupling. In this communica- 
tion, we will use Supernovae data to find confidence re- 
gion for the parameters characterizing the theory. More- 
over, all publications related to chameleon fields used 
quantities in Einstein frame to perform calculations and 
to interpret results. In this communication, we review 



the mechanism by deriving unambiguously physical ob- 
servables from Einstein and Jordan frames. The latter 
is defined as the set of gravitational degrees of freedom 
in which matter couples directly to a metric, although 
that metric does not propagate like it does in general 
relativity. This distinction between Jordan and Einstein 
frame quantities has implications for cosmology because 
we used a Jordan frame conserved densities contrary to 
what is done in [11] [T^. In the same spirit, we take 
a great care in the derivation of the observable (the lu- 
minosity distance). To avoid any confusion, we derive 
observables in both frames showing their physical equiv- 
alence. 

The same approach is applied for computing the 
static and spherical solution representing the solar 
system. Indeed, we use Jordan frame quantity in the 
field equations contrary to what was done in [8l [Ml [15] 
(but following the work of [IS]). On solar system 
scales, analytical solutions representing the scalar field 
were computed in [S] H] [TU [T31 HI] assuming a lot of 
hypothesis that have never been checked. In this paper, 
we compare these solutions with a numerical integration 
of the Einstein equations (in the same way as in 161 ). 
With these solutions, we derive the post-newtonian pa- 
rameters taking care of the gauge used and we compare 
our results with the literature. Finally, we present a 
combined analysis of the cosmological study and the 
post-newtonian constraints. 

The strategy followed is to fit the model on real data 
at two different scales: on cosmological scales and on 
solar system scales. In this paper, we decide to focus 
only on the original theory proposed by Khoury and 
Weltman [5J [3] which is a tensor-scalar theory with an 
exponential coupling function ^(0) — e^'^ and a runaway 
potential V{(t)) of the Ratra-Pcebles type. This model is 
fully described in Section jll} We let for future work the 
study of other coupling constants and/or potentials. 

Concerning cosmological scales presented in Sec- 
tion |III[ we perform a likelihood analysis of the latest 
Supernovae la data. For this analysis, the cosmological 
evolution equations are derived and the expression of 
the luminosity distance is given. In the derivation of 
these equations, the use of the Einstein or Jordan frame 
quantities is fully discussed. In particular, the derivation 
of the luminosity distance has been performed in both 
frames showing their physical equivalence. From the 
likelihood analysis, we find confidence regions for the 
different parameters involved in the model. For models 
within these confidence regions, we also perform a full 
cosmological analysis consisting in a detailed study of 
the evolution of the scalar field 0, the evolution of the 
cosmic expansion and its acceleration and we derive 
an effective dark-energy equation of state from the 
Friedmann-Lemaitre equations. With this analysis, we 
can clearly identify models explaining the accelerated 
expansion of the universe. 
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In Section IV we focus on solar system constraints. 
Solar system constraints are mainly composed of three 
tests: violation of the weak equivalence principle, fifth 
force constraints and Post-Newtonian constraints. In 
our case, the weak equivalence principle is not violated 
since the coupling constant between the scalar field and 
matter is universal. However, we emphasize the fact 
that string theory can give raise to tensor-scalar theory 
with different coupling constants for the different matter 
fields. In this case, violation of the equivalence principle 
is expected [T7]. 

In the model considered here a fifth force term is 
present. A fifth force analysis can be found in Khoury 
and Weltman [8] while we concentrate our attention 
on the analysis of the post-newtonian parameters. To 
derive the expression of the post-newtonian parameters, 
we use the analytical solutions of the field profile in a 
static spherical configuration that can be found in the 
literature [HI HH UHl HB] • First of aU, in order to check the 
consistency of the analytical expressions, we compare 
these analytical solutions (obtained by assuming a lot of 
hypothesis) with a full numerical determination of the 
field profile solving the Einstein equations representing a 
spherical and static body. Then, a careful derivation of 
the post-newtonian parameters is performed. With the 
models within the confidence region of the cosmological 
analysis, we compute the post-newtonian parameters 7 
and we check if they satisfy the constraint ([T]). 

With these two analysis (cosmological and solar system 
analysis), it is possible to answer to the question: can 
the chameleon effect reconcile cosmological observations 
with solar system constraints ? For the model treated 
here, the answer to this question is negative. Indeed, 
we will show that on one hand large coupling constants 
can passed the solar system constraints while they are 
excluded by usual analysis. On the other hand, tensor- 
scalar theories can indeed explain cosmic acceleration. 
But the main conclusion of this communication is that 
unfortunately, the parameters confidence region needed 
to explain cosmic acceleration does not overlap the one 
coming from solar system constraints. This means it is 
not possible to explain cosmic acceleration with the same 
parameters as the one passing solar system experiments. 



II. THE MODEL 

The model studied in this paper is the one proposed 
in the original chameleon paper [8l [9] and is given by the 
following action in the so-called Einstein frame {c— 1): 



Idtt 



where rUp = G~^^^ is the Planck mass, g is the deter- 
minant of the Einstein- frame metric g^,^, R is the Ricci 
scalar and are matter fields. In this frame, the scalar 
field (j) interacts directly with matter particles through a 
conformal coupling characterized by the function A((/)). 
Finally V{(j)) is a potential characterizing the scalar self- 
interaction. This action generalizes the Jordan-Fierz- 
Brans-Dicke theory (recovered in the case of a vanish- 
ing potential and of an exponential coupling function 



A(0) 



) [2111] but also the quintessence model [I91I20] 



which can be recovered in the case A{(j)) = 1. In this com- 
munication, we consider that all matter fields couple in 
the same way to the scalar field with only one coupling 
function A{(j)). The above action can be justified by low- 
energy limit of string theory or supergravity 171 [^THM] . 

Observable quantities are not directly obtained in this 
frame because physical measurements are performed us- 
ing rods and clocks built upon matter fields '^m- These 
matter fields are universally coupled to the so-called Jor- 
dan frame metric 



(3) 



and therefore follows geodesies corresponding to and 
not to the Einstein frame metric 5^1^. One can perform 
the conformal transformation in the action ([2]) to obtain 
the Jordan frame action (see for example [51 lol E51 126| ). 
Nevertheless the Jordan frame is not convenient for 
studying the dynamics of the scalar field and of the met- 
ric 5^1. because the tensor and scalar modes are kinemat- 
ically coupled and because this frame introduces some 
fictitious singularity [25l [26| . 

In the following, a ~ will denote quantities expressed 
in the Jordan frame. Observable quantities can easily 
be computed in this frame in the same way as in Gen- 
eral Relativity (GR). This does not mean that observable 
quantities can only be derived in the Jordan frame. In- 
deed, physical predictions are not dependent of a change 
of variables in terms of a conformal transformation. If 
the derivation of observable quantities in Jordan frame is 
easy, one has to be more careful when working with the 
Einstein frame. In particular, all physical units have to 
be scaled with the coupling function A{ip). 

The potential V{(j)) in the action ([2]) is a key element 
in the chameleon theory. It has to be of the runaway 
form [5J H] . The potential used here is the Ratra-Peebles 
(RP) inverse law potential [HI [20] 



(4) 



where a and A are free parameters. 

In this communication, we only consider an exponen- 
tial coupling function 



A(0) 



(5) 



+5™ [ 



(2) 



This coupling function is the one used in the original 
paper [8, 9 and we let the study of other type of coupling 
function to future work. 
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From variational principles, we can derive field equa- 
tions from action Varying this action with respect to 
the Einstein frame metric gives the Einstein field equa- 
tions 



Rn 



Stt 



-4715^,9,05"^ - -^9^uV{(j}) (6) 



where T^, = -{2/ ^f^){dC^f^ / dg'''') is the Einstein 
frame stress-energy tensor. 

The variation of action ^ with respect to the scalar 
field (t> leads to the Klein-Gordon equation for the scalar 
field: 



fc(0) 



(7) 



where T = g^'^T^^ is the trace of the stress-energy tensor, 
□0 = g^'^V ^4) and k{(f)) is the scalar coupling strength 
to matter and is given by 



d\nA((j)) _ 1 

d4 ~ A{(l)) del) 



(8) 



Finally, the invariance of action ^ under coordinates 
transformations gives the following conservation equation 
in the Einstein frame: 



k{cj))Td''(i). 



(9) 



In the following, we will always consider matter to 
be described by a perfect fluid. Therefore, the Einstein 
frame stress-energy tensor is expressed as 



T^"' = {p + p)u''u'' + pgf" 



(10) 



where p and p are the Einstein frame density and pres- 
sure and is the 4-velocity of the fluid. As can be 
seen from Q, this stress-energy tensor is not conserved 
in the Einstein frame because of the explicit coupling be- 
tween matter and the scalar field. On the other hand, 
the Jordan frame stress-energy tensor is conserved be- 
cause matter field is universally coupled to the metric 
5^11/. Moreover, Jordan frame is defined as the frame in 
which matter experiences locally the laws of special rel- 
ativity. The Jordan frame matter density and pressure 
{p and p) are directly observable. The observable stress- 
energy tensor (Jordan frame) is related to its Einstein 
frame counterpart by a conformal factor [SI [5S] 

T%=A4(</,)f^ , 

= A\cf,){p + p)u^^u, + A*{c^)pS^: 
= AHmp + PWu,^ + A%cj))p6^. 

This implies the following relations between matter den- 
sity and pressure in Einstein and Jordan frame: 



P = A*{c^)p 
p = A^{^)p. 



(11a) 
(lib) 



In the following sections, we will derive observational 
constraints on different parameters characterizing the 
theory (namely parameters involved in the potential and 
the coupling function) and characterizing the cosmologi- 
cal content of the universe. First, we will perform a least- 
square fit on the SuperNovae la data. This involves the 
derivation ant the integration of the equations describ- 
ing the cosmological evolution of the universe. These 
equations will be derived from field equations ([6| and 
([t]). Secondly, we will use solar system observations to 
obtain a different constraint on the theory parameters. 
To do this, we will derive the static spherical solution to 
the field equations. This solution will represent the solar 
system space-time. 



III. COSMOLOGICAL CONSTRAINTS 

This section is devoted to the study of the cosmological 
dynamics of the model described in the previous section. 
Let us briefiy outline our original contributions concern- 
ing chameleon cosmology. We derive the evolution equa- 
tion in Einstein frame but using density and pressure 
defined in Jordan frame. Therefore, we introduce a new 
variable f2 representing the observable density parame- 
ter. In Sec. |III B{ we derive a luminosity distance relation 
which is frame-independent (as any observables should 
be) and therefore allows us to make cosmological pre- 
dictions without ambiguities. In particular, we present 
a likelihood analysis of the latest SNe la data's. These 
analysis extends previous works on quintessence models 
(characterized by a vanishing coupling constant k — 0) 
to the cases of A: 7^ where no analysis exist to our 
knowledge. We identify models statistically consistent 
with data and show show that the statistical adequacy 
decreases with k. The rest of the section comprises a 
detailed analysis of the model dynamics within the confi- 
dence regions previously established. First of all, we ex- 
tend the work of Brax et al [Il| concerning the study of 
the dynamics of the scalar field by introducing three time- 
scales and we showed how the dynamics of depends on 
the parameters of the model. In Sec. |III E we propose an 
original analysis of the cosmic expansion within the mod- 
els explaining data. In this analysis, contributions due to 
the quintessence potential and to non-minimal coupling 
are clearly identified and compared. Finally, we interpret 
considered models in standard FLRW cosmology in Gen- 
eral Relativity with an effective fluid whose quintessence 
and non-minimal coupling contributions to the effective 
equation of state are identified. 



A. Evolution equations 

If we assume a flat Friedmann-Lemaitre-Robertson- 
Walker (FLRW) background space-time, the Einstein 
frame metric can be written as 



-dt^ + a^{t)d£^ = a^{t){-dTf + df) (12) 
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where a{t) is the Einstein frame cosmic scale factor, t is 
Einstein frame cosmic time and r] is conformal time. The 
corresponding Jordan frame metric can be written as 



dS^ = -dp + d^{t)df = A^{<P)ds^ 

= -A^i^{t))dP +A^i^{t))a\t)de^■ 



iU) 



The observable cosmic time t and the observable cosmic 
scale factor (a) are obtained from the last relations: 



di ^ A{(j){t))dt 
a{i) = A{(j){t))a{t)- 



(14a) 
(14b) 



Replacing metric ( 12 ) in Einstein field equation s ([6[ ) and 
using the expression of the stress-energy tensor (111, one 



obtains the Friedmann and acceleration equations 

2 



47r 



'2 



3m2 



aV(0) 



Stt 
Snip 



a'^A^{(f>)p (15a) 



a 
a 



47r 



16n 



Snip 



(15b) 



where a ' denotes a derivative with respect to confor- 
mal time T] and where p and p are the observable matter 
density and pressure. The Klein-Gordon equation ([T]) 
becomes 



-'^a'A\m-m^ 



a" dV 
d(f) 



Finally, the conservation equation ^ gives 



Introducing the equation of state 



P 
P 



P 

~ 7 

P 



(16) 



(17) 



(18) 



and replacing the Einstein frame matter density and pres- 
sure by their Jordan frame counterpart (111, the conser- 



vation equation (17) becomes 

(A(0)a)3(i+-)p : 



Po 



(19) 



where subscript refers to the actual epoch characterized 
by a flo = A{(j)o)ao — 1. From the last equation, we see 
that the Jordan frame stress-energy tensor is conserved 
since we find the usual expression 5,'^''^"'"'^)^ = est. 

We can introduce the observable density parameter de- 
fined by 



n = 



SttGp 



(20) 



where H is the observable Hubble parameter defined by 

(21) 



d dt 



and G is the varying gravitational strength which is re- 
lated to the Einstein frame gravitational constant by p7] 



A\c^) 



mt 



(22) 



Replacing po in the conservation relation ( 19 ) by its 
value derived from ( 20 1 , one finds 



87rA2(0o)(A(0)a)3(i+^)' 



(23) 



We can now use the equation of state ( 18 1 and the ex- 



pression of p (23 1 in the Friedmann, acceleration and 
Klein-Gordon equations ( T5p6 ). If we introduce two 
components in the universe: pressureless matter (rep- 
resenting baryonic and dark matter) characterized by 
u!m = and radiation characterized by = we find 
the cosmological evolution equations 

a 



47r 



12 



Stt 
3m2 



a^V{4>) 



Hi 



^moA{(t)) 



^2(0o) 
47r 



2a 



167r 



aV(0) 



3fc('^) 
8^ A^{d^o) 



rn^ dcf) 



(24a) 



(24b) 



(24c) 



It is worth to comment these equations. First of all, they 
depend on different parameters: the present value of the 
densities parameters ilmo ^nd il^g, the present value of 
the Hubble parameter Hq, the couphng function A{(f)) (in 
this case parametrized by the value of the coupling con- 
stant k) and the potential V{(j)) (in this case parametrized 
by the constants A and a). The quintessence cosmologi- 
cal equations are recovered in the case of a constant cou- 
pling constant {A{(j)) = 1 and k{(j)) = 0) [28]. The above 
equations are somewhat different on what can be found in 
the literature [TTl [^j because we have decided to work 
with Jordan frame density and pressure (these are di- 
rectly observable) following what is done by Damour et 
al. in Uni and in Babichev et al. [inj while most of 
chameleon papers are using some hybrid conserved den- 
sity [p) that can be related to the traditional Einstein 
frame density [p non-conserved due to the coupling to 
the scalar field, see ([9])) or to the Jordan frame density 
pressure (p conserved in Jordan frame) [5J IHl HH US] : 



p = A-\ci>)p = A^{^yp- 



(25) 
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Naively, all these definitions of matter density could be 
used in the source terms of the field equations. Never- 
theless from an interpretative point of view, the use of 
the Jordan frame density p is justified since Jordan frame 
quantities are measured as in GR with metric 5^1/ cou- 
pling universally to matter fields. Moreover, equation 
(251 is not robust since it depends on symmetries of the 



cosmological principle and, even in that case, is only valid 
for pressureless matter. Jordan Frame quantities are con- 
served ones, while not suffering from these pathologies, 
and can further safely be used in solar system physics. 



B. Luminosity distance 

SNe la standard candles provide a measurement of the 
observable luminosity distance ^^(z). This luminosity 
distance is easily obtained in the Jordan frame. Never- 
theless, to illustrate the fact that observable quantities 
can also be derived in Einstein frame, we will derive the 
expression for in both frames. The present result is 
therefore frame independent. 

In Jordan frame, the luminosity distance is expressed 
as in GR 



dL{z) = {l + z) 



dy 

H{y) 



(26) 



where z = l/a--lis the observable cosmological redshift. 
Finally, the distance modulus is defined as 



/i = 25 + 5 log 



dL 



1 Mpc 



(27) 



If we want to express the same observable quantity from 
Einstein frame, we introduce the Einstein frame distance 
luminosity 



dL{z) = il + z) 



dy 
H{y) 



(28) 



where z — 1/a — I and H is the Einstein frame Hubble 
parameter H — = In Einstein frame, the physi- 
cal units are rescaled by a factor A^cj)). For this reason, 
the observable distance modulus is given by 



^ = 25 + 5 log 



dL 



A{4)) X 1 Mpc 



(29) 



The equivalence between the two expressions (27) and 



(291 can be shown. Substituting I + z — 1/a, H 
and dz = in we find 



dL{z) 



and thus 



1 f"" da 1 I ' . 770 — 77 



a^H 



^ = 25 + 5 log 



drj 



Vo 



(30) 



(31) 



(if ?7 is expressed in MPc which is possible since C77 has 
dimension of a length). 

On the other hand, using relation 
that dt — adrj we can express the Hubt 



14a) and the fact 



Die parameter as 



a dt 



1 dd 



dA((j)) dt d 



n2 



Substituting this relation and using the fact that dz 
— , the relation ( 26 ) becomes 



MS) 



dd 



drj 



Vo - V 



(32) 



Vo 



Finally, the distance modulus (27) is unambiguously 
given by 



/2 = 25 -f- 5 log 



25 + 5 log 



Vo - V 
A(0)a 



with 7] expressed in MPc. Relation (Mb]) shows that 
the last expression is equivalent to (31). This shows the 



physical equivalence between the two frames and gives a 
useful formula to compute the distance modulus without 
any ambiguities from any frame. The distance modulus 
can be evaluated from the integration of the cosmological 
evolution equations ( 24 ) that gives the evolution of a and 
(f) with respect to rj. 



C. Supernova likelihood analysis 

In this section, we present the results of a likelihood 
analysis of recent SN la data [30] . Our goal is to identify 
models characterized by a set of parameters which are 
statistically consistent with observations, while departing 
from the standard ACDM values. Cosmological models 
are characterized by 6 parameters: flmo, ^ro, Ho, k, A 
and a. Since the present radiation density parameter is 
very low, it has a negligible influence on SN la measure- 
ments. Therefore, we fix its value Ctro — 7.97 10~^ [30j . 
The value of the energy scale of the potential A is opti- 
mized such that for a given value of a, the input value 
of flrno is retrieved. Finally, since the Hubble diagram 
leaves Hq unconstrained in the plane a — rimo |3H 132) , 
we fix the value of Hq ~ 70 A km/ s /MPc [33]. Because of 
the degeneracy of the Hubble diagram with respect to Hq, 
we compute the marginalized likelihood over Hq |34l 135) . 

After the above considerations, each model is finally 
characterized by a set of three parameters {i~lmO, ct, k). 
We run a series of model decomposing the parameters 
domain using a uniform grid. The results of this analysis 
is summarized on Fig [T] where we plot the conditional 68 
% and 95 % confidence regions in the plane ^moh^ ~ 01 
(with h = km/s/Mpc) with an assumed value of 

the coupling constant k. As can be seen, the concordance 
ACDM model (corresponding to /c = and a = 0) is 
within this confidence region. Models characterized by 
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0.15 



0.15 




0.15 



mO 



mO 



FIG. 1. Representation of the conditional 68 % (in blue-dark) and 95% (in green) confidence regions in the h^tlmo — ol plane 
for different assumed values of coupling constant k from the analysis of the UNION SN la Hubble diagram. The star represents 
the best fit. 



k — Q are quintessence model. As noticed from previous 
quintessence data analysis |281 136j . quintessence models 
tend to fit the data by requiring values of flmoh^ lower 
than in ACDM. This trend is reinforced when increasing 
the value of the coupling constant k. In addition, we can 
see that the area of the confidence region does not reduce 
significantly with an increase of the coupling constant k. 
This means that the fine-tuning on parameters a and 
flmO does not evolve with k. This is mainly due to new 
dynamical regimes explaining cosmic expansion with a 
small value of the matter density parameter but a high 
value of the coupling constant as can be seen further. 

Fig. [2] shows the evolution of the with the observed 
cosmological density fi„io for a fixed value of a = 0.5. 
The gives an idea of the goodness of the fit of the 
model to the data. An increase of the coupling constant 
k results in an increase of the minimum value of the 
which shows that models with non-vanishing coupling 
constant provide a slightly worst reproduction of data. 
For high values of the coupling constant, the curve of the 
has two minima: one with a small value of fimo and 
one with a usual value of r2„iO- These two minima char- 
acterize two very different cosmological evolutions (as it 



will be discussed below) and they lead to very different 
prediction for the age of the universe. The prediction for 
the age of the universe is represented on top of Fig. [2] as a 
function of O^g. First of all, it is interesting to note that 
models with a high coupling constant lead to a higher es- 
timation of the age of universe than quintessence models 
with the same cosmological parameters. For example for 
h^ttjno = 0.1, the estimation of the age of the universe of 
a quintessence model is 14.2 10^ years while for the same 
parameters, the model with a coupling constant fc = 10 
gives an age of 15.3 10® years. Moreover, models with a 
small value of (imo (allowed by the likelihood analysis) 
leads to an estimation of the age of the universe much 
higher. 

A supplementary interesting result is the relationship 
between A and the other parameters. Fig. [3] represents 
the evolution of A with respect to a for different values 
of k and ^mo- It can be seen that the influence of the 
coupling constant and of the matter density parameter is 
very weak in comparison with the influence of a. We flt 
a curve on these data and we flnd the following relation 



log A 



19a - 47 



(33) 
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FIG. 2. Top; Representation of the age of the universe with 
respect to the observed matter density ClmO for models char- 
acterized by Q = 0.5. 

Bottom: Representation of the chi squared per degrees of free- 
dom judoS (number of degree of freedom) with respect to 
the observed matter density f2mO for models characterized by 
a = 0.5. 



which is the same relation as found in |32] for 
quintessence model. The last relation gives an order of 
magnitude for the value of A. 



D. Evolution of the scalar field 

It is interesting to study the cosmological evolution of 
the scalar field for cosmological model compatible with 
the SN le analysis from previous section. The Klein Gor- 
don equation governing the scalar field evolution is given 
by ( 24c ) and can be written as 



2li(p = -fc(0) — p„ 



cff 



d4> 



(34) 



where Ti. — a' /a and T4ft(0) is an effective potential 
whose expression is 



1 



(35) 



With the hypothesis adopted in this paper (the use of 
a Ratra-Peebles potential and an exponential coupling 
function), this effective potential exhibits a minimum 
(see Fig. 10). The scalar field is thus attracted toward 
the minimum of this effective potential. This minimum 




FIG. 3. In blue: data representing A with respect to a for 
different value of k and Q.mO- It can be seen that the influence 
of k and Clmo is very weak compared to a. In red: curve 
resulting from a fit on these data (see (33 1). 



depends on the scale factor a (due to the fact that the ef- 
fective potential depends on p^) and therefore is moving 
continuously with time. The shape of the effective poten- 



tial for two different densities is represented in Fig. 10 
The density is decreasing with the scale factor (see equa- 
tion ( |19[ )) which implies this effective potential was very 
narrow in the early time and it becomes wider and wider 
with the time. 



It is worth noticing that the expression ( 35 1 of the 
effective potential is somewhat different of the expression 
given in Brax et al |llj . The difference comes only from 
the fact that we use the Jordan frame density p while 
Brax et al used a hybrid definition of the density p given 
by (|25|. 

Working with Ratra-Peebles potential V{4)) = 
lm'^(j)°^ and with an exponential coupling function 
A{(j>) = e^"^ gives an analytical expression for the mini- 
mum of the effective potential determined by the follow- 
ing conditions 



= 



vA4+" 



kpr, 



(36) 



Solving this equation gives the value of the scalar field 
that minimizes the effective potential (pmm as a function 
of the cosmological density Pm, of the coupling constant 
k and of the parameters of the potential {a and A): 



with b = 
tion [37J. 



{Prn) 



1 



4fc 



LW 



4kb 



(37) 



and LW{x) the W-Lambert func- 



9 



For models fitting the SNe la (see previous section), 
the typical behavior of the scalar field if the following: 

• at the beginning, the attractor mechanism is very 
efficient and the scalar field is attracted by the 
minimum of the effective potential and oscillates 
around it. This is due to the fact that the effective 
potential is very narrow in the early time. Fig. [4] 
shows the behavior of the scalar field for different 
initial conditions and the minimum of the effective 
potential (pmin (37). We can clearly see that the 



attractor mechanism is very strong for small scale 
factor. This implies that the cosmological evolu- 
tion is very weakly sensitive to the scalar field ini- 
tial conditions (see Fig. |4]) . This is a reminiscence 
of the tracking properties of the Ratra- Peebles po- 
tential [3S]. 

• during some time, the scalar field follows closely 
the minimum of the effective potential. 

• depending on the model considered, the attrac- 
tor mechanism may be not strong enough and the 
scalar field may not be able to follow the minimum 
of the potential (see Fig. [4] and Fig. |6|. In par- 
ticular for low coupling constants, the effective po- 
tential becomes too wide to attract sufficiently the 
scalar field. 

The necessary condition for the scalar field to follow 
the minimum of the effective potential can be derived 
by considering the time scales involved in the scalar field 
evolution. There are three different time scales: 

• the damping time-scale due to the in equation 
( 34 1 term given by 



1 



2n 



= 2H 



(38) 



amp 



where tdamp is an Einstein frame cosmic time-scale 
characterizing the Hubble damping. 

• the time-scale characterizing the scalar field oscil- 
lations around the minimum of the potential. This 
time scale is related to the effective mass of the 
scalar field given by 



m 



off 



The related time-scale is given by 



1 



(2^) 



(39) 



(40) 



which characterizes the response time of the scalar 
field. 

the last time-scale involved in the scalar field evo- 
lution is related to the evolution of the effective 
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FIG. 4. Evolution of the scalar field as a function of the 
Jordan frame scale factor a (the observable one) for differ- 
ent initial conditions ij>i and for the model characterized by 
(fe = 1, fimo = 0.197 and a = 0.5). The minimum of the 
effective potential (see relation ( |37[ )) is represented in blue. 
On top: curves on a log-log scale, we can see that the scalar 
field reaches the minimum of the effective potential for all the 
initial conditions. On bottom: linear x-axis, we see that the 
scalar field does not follow the minimum of the effective po- 
tential from a ~ 0.1 and the evolution of the scalar field is 
independent of the initial condition (f>i (all curves are super- 
imposed). 



potential. This potential is continuously moving in 
time due to the presence of the pm term in ( 35 1 . 



We characterize the time-scale due to this evolu- 
tion by the time scale related to the variation of 
the minimum of the potential 4>min 



1 



1 



dp, 



da 



(41) 



where a dot denotes the derivative with respect to 
Einstein frame cosmic time [dt). Using the expres- 



sion of I 



a ( 14b ) , we finally get 

1 sfi? 



(37), of Pm (23) and the definition of 



1 



4fc(/)„ 



(42) 



The scalar field follows the minimum of the effective 
potential if the damping time is higher than the scalar 
field response time {tdamp > t^) and if the evolution of 
the effective potential is slower than the field response 
time {tmin > ttj))- In other words, the conditions charac- 
terizing the fact that the scalar field follows the minimum 



10 



of the potential are given by 



> 



m 

4^?2 > (T 



a 



(43a) 
(43b) 



The three time-scales are represented on Fig. [5] (for a 
model fitting SN la data with fc = 1). We can see that the 
intersection of with the other curves occurred at a ^ 
0.1 which is exactly the scale factor where the scalar field 
stops to follow the minimum of the effective potential (sec 
Fig. [4]). We note at the Jordan frame scale factor from 
which the scalar field stops to follow the minimum of VoS- 
Fig. [6] (a) represents the evolution of the scalar field and 
of the minimum of the effective potential for different 
values of the coupling constant k. The transition scale 
factor (it increases with k and for k ^ 15, the scalar 
field stays in the minimum of the potential during all the 
evolution. In Fig. [6] (b), the behavior of the scalar field 
for the second minimum of the curves for fc = 10 is also 
represented. The evolution of this model (characterized 
by a small value of the density parameter r2„io) is very 
different from the other ones. The scalar field increases 
very fast. 

Fig. [6] (c) represents also the evolution of the transition 
scale factor with the coupling constant fc {at being com- 
puted by the comparison of the Hubble damping time 
scale and the field response time). The transition scale 
factor increases with the coupling constant which is log- 
ical since the effective potential is tighter for high cou- 
pling constant, it is therefore more difficult to leave this 
potential. 
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FIG. 6. (a) : Evolution of the scalar field (continuous lines) 
and of the minimum of the efi'ective potential (j)min (371 
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(dashed lines) for diflferent values of the coupling constant 
fc. For small values of the coupling constant, the scalar field 
leaves 4'min quite early. The transition scale factor at where 
the field leaves the potential increases with the coupling con- 
stant fc. For fc = 15, the scalar field is still in the minimum of 
the effective potential at the present epoch. 

(b) : Evolution of the scalar field (continuous lines) and of the 
minimum of the effective potential 4>m.in ( |37| l (dashed lines) 
for two models with fc = 10 located in the local minima of 
the (see Fig. [2]). The evolution of the model located in the 
second minimum of the located in the low fimo minimum 
is very fast. 

(c) : Evolution of the transition scale factor at with the cou- 
pling constant fc. 



FIG. 5. Evolution o f th e three time scales chara cter izing the 
field response time (40 1, the Hubble damping (38 1 and the 



evolution of the effective potential (41 1 as a function of the 



Jordan frame cosmic scale factor 5 for the same model as in 
Fig. |4] We can see the intersection between and tdamp 
occurred at the transition scale factor fit that corresponds to 
the scale factor where the field leaves the minimum of the 
potential (see Fig. [4|. 



E. Evolution of the cosmic expansion 



We study here what is the origin of the observed cos- 
mic acceleration. In this context, cosmic expansion is de- 
scribed in the Jordan frame by the scale factor a = A{<j))a 
that is measured with matter. In order to identify the 
origin of the acceleration, we need t o com pute the value 
of 



1 cPa 
a dt^ 



With the definition of a ( 14b ) and with the 
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evolution equations ( 24 1 , we can express 



1 d^h 



where all the terms are given by 
2 



Qm + Qq + Qnmc 



Qm = 



1Q 



Stt 



(44) 

(45a) 
(45b) 



9NMC 



(45c) 



3m2A2(0)i72 
k(l>" 

The first term gives the usual term present in GR involv- 
ing matter (here the subscript M refers to presure-less 
baryonic and dark matter and also to radiation). The 
second term qq is a term due to the presence of the 
scalar field in the dynamics. This term is also found 
in quintessence models (where the coupling function is 
equal to A(^) = 1 and the coupling constant k van- 
ishes) [32l |35j EB]. Finally, the last term q-NMC is due 
to the non-minimal coupling (A(</)) ^ est) and is closely 
related to the evolution of the scalar field. 

The different contributions to the cosmic acceleration 
are represented on Fig. [7] for three different values of 
the coupling constant k {k — 1, k — 10 and k — 20) 
for best-fit model. The matter contribution is more or 
less the same for the different k. What is interesting 
to notice is that for small value of fc, the cosmic accel- 
eration is explained by the quintessence term qg while 
the non-minimal coupling does not play any role. An 
increase of the coupling constant produces a decrease of 
the quintessence contribution in favor of the non-minimal 
coupling contribution ^nmc- In conclusion, if k » 1, 
the cosmic acceleration is explained by the non-minimal 
coupling. 

In order to compare the difference in the dynamic of 
the two minima in the curves (see Fig. [2] (bottom) and 
discussion related). Fig. [t] represents also the evolution of 
the acceleration factors for the two different minimum for 
k = 10. The evolution of the model in the left-minimum 
(with small matter density parameter) is very fast and 
reaches an asymptotic behavior with the quintessence 
contribution being of the same order of magnitude than 
the contribution from the non minimal coupling. The 
matter contribution tends to a very small value which is 
logical since the density parameter of this model is very 
small. 



F. Dark-Energy effective equation of state 

Following what is usually done in quintessence stud- 
ies m [21 133 [Ml and in the context of the AWE 
hypothesis [40] , we can modeled the presence of the non- 
minimally coupled scalar field by an effective fluid of 
dark energy in GR (parametrized by a density param- 
eter fluE and by the equation of state parameter ujde)- 
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FIG. 7. Representation of the different terms of the acceler- 
ation factor for different coupling constant fc. The different 
terms qm for the matter contribution, qq for the quintessence 
contribution and ^nmc for the non minimal coupling contri- 
bution are given by equations ( 45 I . The three first curves 
fc = 1, 10, 20 represents model in the global minimum of the 
while the last curve k = 10— left represents a model in the 
second minimum (left) of the (see Fig. [2]). For model in the 
global minimum of the x^- the GR contribution is of the same 
order of magnitude for different values of k; the quintessence 
potential contribution decreases when the coupling constant 
increases to leave the room for the contribution coming from 
the non-minimal coupling. The model in the low ClmO min- 
imum of (left) characterized by a lower value of the QmO 
has a much faster evolution. 



In order to identify clearly the quintessence contribution 
from the contribution coming from the non minimal cou- 
pling, we introduce two effective fluids in GR related to 
quintessence term (characterized by flq and ujq) and to 
the non minimal coupling (characterized by JInmc and 

'^NMc)- 

The observed Hubble constant is defined by (21 1. Us- 



ing relations (14b I and (14a I, we can relate this Jordan 



frame Hubble constant to the Einstein frame one 



H = 



1 da 
a di 



1 d{Aa) 
Aa Adt 



1 



H + a((l)) 



d<j) 
'dt 



(46) 



Inserting the evolution equations ( 24 1 in ( 46 1 and by 
matching it with the usual Friedmann equations in the 
Jordan frame 

1 = n„, + f^r + ^Q + f^NMC, (47) 

one identifies the expression of the observable density pa- 



12 



rameters 



Stt 



The total parameter density is given by 



1/(0) (48a) 



(48b) 



DE 



fir 



n 



NMC 



(49) 



Fig. |8] (a-b) represents the evolution of the different den- 
sity parameters fl defined above for different coupling 
constants k. The quintessence density parameter de- 
creases when k increases while the non-minimal coupling 
density parameter increases with k. This is consistent 
with the behavior of the acceleration factor q illustrated 
in Fig.[7j For high value of the coupling constant, the cos- 
mic acceleration is mainly explained by the non-minimal 
coupling. The asymptotic behavior of the left-minimum 
model (low value of fimo) is also confirmed by Fig. [8](b). 



Matching the acceleration equation ( 44 ) with the usual 
GR acceleration 

q = ^ — f2r " 2^'^'-"'^ ^ 3a)g) — -f^NMc(l + Scjnmc) 

gives the equations of state 

m2072 - V{(j>) 



'^'^ m2072 + 1/(0) 

3 2^+ k4>^ 
and the total equation of state is given by 



wnmc 



(50a) 
(50b) 



^DE 



(51) 



DE 



Fig. |8] (c) represents the evolution of the total equation 
of state parameter Qde for different coupling constants 
k. All models represented on this figure are statisti- 
cally equivalent to ACDM for the SNe la measurements. 
An increase of the coupling constant produce a signi- 
ficative change in the shape of the curves. In particu- 
lar, the value of can become positive (this deriva- 
tive is sometimes noted —uji in a Chevallier-Polarski- 
Linder parametrization of the equation of state param- 
eter uj{a) = 0)0+0)1(1 — 0,) [U 112]). In quintessence 
scenario with Ratra-Peebles potential, the derivative of 
0) is always negative (wi > 0) [JS]. Therefore, a mea- 
surement of the value of the derivative of the equation of 
state parameter should in principle discriminate strongly 
non minimally coupled DE and uncoupled quintessence. 
Such a measurement will be performed by the EUCLID 
mission recently selected by ESA. 
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FIG. 8. (a-b) : Evolution of the density parameter fl corre- 



sponding to the quintessence contribution {0.q given by ( 50 
and represented in blue (dark)) and to th e non- minimal cou 
pling contribution (S^nmc given by j50b| and represented in 
red (light)) different coupling constant: Fig. (a) : k = 1 and 
fc = 20 ; Fig. (b): two models characterized with fe = 10 and 
in the different local minima of the 

(c) : evolution of the equation of state parameter lode as 
a function of the cosmic scale factor for best-fit model with 
different coupling constants and for a model in the second 
minimum of the with k — 10. 



IV. SOLAR SYSTEM CONSTRAINTS 



In this part, we will study the constraints on the pa- 
rameters characterizing self interaction of the scalar field 
and on the coupling between matter and the scalar field 
that can be obtained from solar system experiments. Tra- 
ditional tensor-scalar (traditional in the sense where no 
potential is considered or chameleon mechanism is not 
playing any role) is severely constrained by solar sys- 
tem and by binary pulsar experiments. For example, 
deviations from General Relativity can be expressed in 
the Paramctrized-Post-Newtonian (PPN) framework [1]. 
The PPN parameters related to a traditional tensor- 
scalar theory of gravity with an exponential coupling con- 
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stant A^cf)) = e'"^ are given by (see e.g. [S]|^ 



7-1 



/3 - 1 = 0- 



(52) 



The present constraint on the 7 parameter is obtained 
from the Shapiro measurement of the Cassini probe 



7 - 1 = (2.1 ±2.3) X 10~ 



(53) 



This constraint on 7 shows that tensor-scalar theory of 
gravity can be viable only for very small linear coupling 
parametrized by k. The main point of the chameleon 
fields presented in Khoury and Wcltman [8j [9] is that 
this constraint can be evaded thanks to the chameleon 
mechanism. If the potential of the scalar field is chosen 
so that its effective mass becomes large in presence 
of matter, the coupling constant k is replaced by an 
effective coupling constant fceff that can be strongly 
reduced with respect to k. The goal of this part of the 
paper is to study a simplified model of Sun and to derive 
constraints so that the chameleon mechanism is strong 
enough to pass the PPN test of gravity. Furthermore, we 
will compare the admissible parameters region with the 
confidence region obtained by a cosmological analysis in 
the previous section. 

Let us summarize the innovative points of this publica- 
tion for what concerns solar system physics. In the next 
section, we study properties related to spherical static so- 
lutions of field equations (6][7|9). The field equations de- 
riving from a Schwarzschild-like metric is derived. These 
equations are similar that the ones found in [T6] but the 
problem of the boundary conditions is fully discussed in 
the light of the cosmological analysis from the previous 
section. In Sec. IV B we review the chameleon mecha- 
nism presented in 8, 9J. Sec. IV C| is devoted to the anal- 
ysis of the analytical field profile derived in the literature 
[SI [m 1121 [H] ■ In particular, the hypothesis done to ob- 
tain these analytical results are discussed and a compari- 
son between the analytical results and the full numerical 
resolution of the field equations is performed in order to 
identify the validity regime of the analytical solution. 

In Sec. |IVDl we derive the post-newtonian parameters 
related to the solar system metric. Different PPN anal- 
ysis of the chameleon mechanism can be found in the 
literature [3 [Ml . Unfortunately, the present results 
are derived using different gauge which can be mislead- 
ing. Therefore, we present a complete derivation of the 
Post-Newtonian Parameters following the procedure de- 
scribed in Damour en Esposito-Farese [5] and we compare 



^ Note that our scalar field <f> is related to the scalar field ip used by 
Damour and Esposito-Farese by the relation = —?=■ There- 

V 47r 

fore, our coupling constant k is related to their constant a by 
k = y/4iTa and the 7 parameter given by 7 = becomes 



it to current results. Finally, we combine the cosmolog- 
ical analysis and the Post-Newtonian analysis to derive 
constraints on the parameters of the theory. In partic- 
ular, we show that the considered model cannot explain 
cosmic expansion and satisfy solar-system constraints at 
the same time. 



A. Static spherical configuration 

The metric characterizing a static spherically symmet- 
ric space-time can be written in Schwarzschild coordi- 
nates (in Einstein frame) 



(54) 



Replacing this expression of the metric in the field equa- 
tions and using the definition of the stress-energy tensor 
( 11 1, we find the following equations 



r 



Stt 



Stt 



re^V{cj)) 



(55a) 

(55b) 
(55c) 



- + -(//' - A') ^ + — 



r 2 
e 

'p 



km~p-ip) 



ml d(j) 



p' - -{p + p) + 



(55d) 
(55e) 



where a prime denotes here the derivative with respect to 
the radial coordinate r. These equations extend the ones 
derived by Damour and Esposito-Farese [TUI (where no 
potential was considered) and are equivalent to the ones 
in Babichev and Langlois [35]. This set of 5 equations 
with 6 unknowns needs to be completed by an equation of 
state. In this communication, we consider a simple model 
of star/planet assuming the energy density is constant 
inside the body (/?&). Outside the body, we consider a 
cosmological background of baryonic gas and dark matter 
whose constant density poo is related to the cosmological 
density parameter flmo by the relation (23). 

To understand the chameleon effect, it is useful to con- 



sider the equations ( 55d I in the non-relativistic limit (ig- 



noring the backreaction of the scalar field and considering 
e'^ ^ 1) and neglecting the pressure with respect to the 
energy density: 
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where the effective potential is exactly the one appear- 
ing in the cosmological evolution ( |35[ ). With a runaway 
potential, this effective potential has a minimum given 
by (t>,nin{p)- In the case of a Ratra-Peebles potential, 
this minimum is given by the expression ( 37 1 . As al- 
ready stated in Section |III the expression of the ef- 
fective potential is different from what is usually used 
in the chameleon literature 2, 8, 9] [TTJ [TH [THj because 
we use the always conserved Jordan frame density p that 
are observable (we follow what is done by Damour et 
al. [HI Uni 123 and by Babichev and Langlois [TB|) while 
most chameleon papers are using an hybrid density for 
pressureless matter that is conserved in cosmological con- 
text (p) related to the Jordan frame or to the Einstein 
frame density by a conformal factor as indicated by the 
relation (25 1. This difference does not change qualita- 



tively the results obtained. 

Five boundary conditions are needed to solve the prob- 
lem given by equations (551 and the equation of state. 



Regularity conditions at the origin of coordinates r = 
impose that 



A(r = 0) = 
-^{r = 0) = 0- 



(57a) 
(57b) 



The matter pressure needs to vanish at the boundary of 
the body considered r — R^, giving the supplementary 
condition: 



p{r = Rb) = 0- 



(57c) 



An additional condition on v is still needed, which cor- 
responds to fixing the time coordinate. This choice can 
be arbitrary and does not change the result of the inte- 
gration since the equations depend only on ly' . 

The last condition concerns the scalar cj). Following 
what is done in Khoury and Weltman [5J IH] , in Babichev 
and Langlois j45j . in Tamaki and Tsujikawa [14 and 
in Tsujikawa et al [151 , we suppose that the scalar field 
reaches the minimum of the potential at very large dis- 
tance of the body: 



(f>{r 



(57d) 



The validity of the assumption is questionable. Indeed, 
we have seen in section [lII D| that the cosmological evolu- 
tion does not always produce a situation where the scalar 
field stays in the minimum of the effective potential. In 
particular, for low values of the coupling constant (fc be- 
low ^ 15), we have shown that this is not the case and 
the cosmological value of the scalar field at the present 
epoch (j)Q is different from the value minimizing the ef- 
fective potential {4>min) as illustrated in Fig. [6] In Fig.jo] 
we represent the ratio (/>oo/0o (with being the value 
of the scalar field minimizing the effective potential). It 
can be seen that for low value of the coupling constant, 
this ratio is different from 1 and the boundary conditions 
(57d) can not be fuUfiUed. To be completely consistent. 



one should solve the cosmological evolution of the struc- 
tures in order to have an idea of the field boundary con- 



ditions at the limit of the solar system. This work is out 
of the scope of this paper. 




FIG. 9. Ratio cf> min{p oo) / 4>o between the value of 4> minimiz- 
ing the cosmological effective potential Ves and 4>q the value 
of the field given by the cosmological evolution. The different 
lines correspond to different values of the coupling constant k 
and the filled areas represent different values of Q.mO- 



The five equations ( 55 ) with the considered equation 



of state and with the five boundary conditions (57 1 de- 
scribed above form a Boundary Value Problem (BVP). 
This BVP has been solved numerically by using multigrid 
methods. 



B. The chameleon mechanism 

As already stated in the previous section, in the non- 
relativistic limit, the evolution of (j) is governed by the 
effective potential Vcs{<t>Tp) which depends explicitly on 
the local matter density. This potential is illustrated on 
Fig. [To] where the blue line represent the effective po- 
tential outside the body (the cosmological density p^o is 
much smaller than the body density ph) and the red line 
represents the potential inside a body (star or planet). 
The minimum of the potential is noted (pmin and the 
mass of small fiuctuations around 4imin is obtained by 
evaluating the second derivative of the potential 



902 



4^4 



)k-p 



(58) 



As can be seen on Fig. 10 larger value of p corresponds 
to smaller value of 4>min and to larger value of m (the po- 
tential becomes more narrow). In other word, the denser 
the environment, the more massive the chameleon. In 
what follow, indices h refers to quantity evaluated at the 
minimum of the effective potential inside the body and 
oo refers to quantity evaluated at the minimum of the 
effective potential outside the body. 
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Low Density 

High Density 




FIG. 10. Representation of the effective potential for large 
and small value of the matter density p. This illustrates that, 
as p decreases, the effective potential becomes wider and the 
minimum shifts to higher value. This situation can represent 
the effective potential inside a body (star or planet) in red 
and outside the body in blue. 



In chameleon model, it is assumed that the scalar field 
reaches its minimum far away from the body {(j){r = 
oo) = 4>ao) as we discussed around (57d). On the other 
side, the field does not necessary reach its minimum in- 
side the body. Two regimes are possible: 



• the field does effectively reach the minimum of the 
potential inside the body {(j){r = 0) w and is 
frozen in this minimum in a large part of the body. 
This regime is called the thin-shell regime [8l|9] be- 
cause the field evolves only in a small shell around 
the body's surface. If we note by Rb the body ra- 
dius and by Rr the radius where the field begins to 
move, the condition to have a thin-shell regime is 
given by 

ARh — Rr 

(59) 



Rh 



Rb 



< 1 



Following Khoury and Weltman [51 [3] , we can write 
this condition as follow 



ARb 

Rb 



(t>b 



47T 



< 1 



(60) 



where e is called the thin-shell parameter and $f, 
is the Newtonian potential at the body's surface. 



^ The definition of the thin-shell parameter e is different from the 
definition given in Khoury and Weltman 0[9] because the start- 
ing action is not the same. In particular, our scalar field <j) is re- 
lated to their scalar field 4>KW by the relation if> = "^^"^ and as a 
consequence, our scalar coupling k is related to their scalar cou- 
pling /3 by fe = \/8tt/3. With these relations, the thin-shell param- 
eters defined in Khoury and Weltman [8] |9] £ = '^-^'g^ff/tf""'' 



becomes 



• the field does not reach the minimum of the po- 
tential inside the body and the field is not frozen 
at all inside the body. This regime is called the 
thick-shell regime. 

Fig. [TT] shows qualitatively the difference between the 
two regimes. In the thin-shell regime, the field is frozen 
around 0b inside the body while this minimum of the 
effective potential is not reached in the thick-shell case. 
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FIG. 11. Representation of the evolution of the scalar field 
as a function of the radius coordinate in the two regimes. In 
red (ligth): in the thin-shell regime, we see that the field is 
frozen in the minimum of the potential inside the star. In 
blue (dark): the thick-shell regime, the field does not reach 
the minimum of the effective potential inside the star. 



C. Field profile and fine-tuning 

The field profile 0(r) has been derived analytically by 
Khoury and Weltman [5]. A more precise solution has 
been derived by Waterhouse (TH] and by Tamaki and 
Tsujikawa [14j . They introduced a third regime: the no- 
shell regime which essentially corresponds to the thin- 
shell regime in the limit where the thin-shell vanishes. 
These solutions have been derived by considering a lot of 
hypothesis: 

• the solution has been derived using a Minkowski 
background space-time. 

• outside the body, the derivative of the effective po- 
tential is linearized 

• in the thin-shell regime, inside the body the deriva- 
tive of the effective potential is also linearized 

• rriaoRb « 1 

• the factor A'*((/)) is always approximated to 1 
(which means that kcj) << 1) 
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The validity of the previous hypothesis have never been 
checked afterwards by comparing the analytical solution 
with the numerical solution of the full relativistic system 
of equations ([55]). Finally, Tsujikawa et al. [TS] derived 
the field profile considering a perturbation coming from 
the gravitational background. 

The exterior solution derived with the previous hy- 
pothesis is given by 



kcs GMb 



An rc? 

where fceff depends on the thin-shell parameter e 



(61) 



(62) 



Two regimes are identified depending on the value of e: 



• if e < i 



This case corresponds to the 



thin-shell regime and the effective coupling con- 
stant is given by 



'Ceff 



3/fc 
( 



Rb 



Rb 



1 

rubRb 

2 Ai?f, 1 
jUbRb Rb [mbRbY 



(63a) 



if e 



This case corresponds to the 



thick-shell regime and the effective coupling con- 
stant is given by 



k- 



(63b) 



The analytical expression for the interior solution can be 
found for example in Tamaki and Tsujikawa [14] but will 
not be needed in our case since we will be only interested 
in the exterior solution (which represents the solution 
around the Sun). 

Fig. [12] represents a comparison between the analyti- 
cal solution and the numerical solution of the full BVP 



problem ( 55 1 for a situation corresponding to the sun and 
with a cosmological density of lQ~^''kg/w? for a model 
giving a thick-shell. We can see that the analytical solu- 
tion is very good far from the sun but this agreement is 
less good inside the sun. This is mainly due to the fact 
that the analytical solution is ob tained by neglecting the 
A^{(j)) term in the equation (55d|. Since, for the determi- 



nation of the PPN parameters, we only need the outside 



solution, the analysis of Fig. 12 shows that the analytical 
solution can be used in the case of thick-shell. 

The case of a thin-shell profile around the Sun is more 
delicate to treat from a numerical point of view. Indeed, 
it is no longer possible to solve numerically the full BVP. 
This is mainly due to the fact that pb/poo ~ 10'^° which 
means that, in order to have a model exhibiting a thin- 
shell profile, the effective mass of the scalar field inside 
the body (m^) needs to be huge (as can be inferred from 
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FIG. 12. Representation of the evolution of the scalar field 
around the Sun for a model characterized by the parameters: 
A = 10"^ GeV, a = 0.5, fc = 1 and = IQ-^^itg/m^ With 
these parameters, we obtain a thick-shell around the sun. On 
top: comparison of the field profile around the neighborhood 
of the Sun obtained numerically by solving the full BVP prob- 
lem ( 55 I and by using the analytical solution. On bottom: 



relative difference between the numerical solution and the 
analytical approximation. One can see that the analytical 
solution is very good far from the sun but is less efficient in 
the Sun. 



the relation (83) of Tsujikawa et al. P^S]). As a result 
of this huge effective mass m;,, variations of the scalar 
field inside the body are extremely small and very often 
smaller than the epsilon machine. This leads to prob- 
lems very delicate to treat numerically. The thin-shell 
profile obtained by solving the full BVP was obtained 
by Babichev and Langlois [IS] but in the case where the 
ratio of density was of the order of 10~^ — 10~^. 

Instead of solving the full BVP, which is very hard to 
do, we treat a simplified BVP constituted of the scalar 
field only. This means that we approximate the space- 
time background by a Minkowski space-time (this is also 
one of the hypothesis done in order to get the analytical 
solution). In this case, the one second order differential 
equations with boundary condition {4>' {r = 0) = and 
(/)(r = oo) = 4>oa) is solved by a shooting method im- 
plemented with a stiff integrator in quadruple precision. 
With this method, it is possible to get thin-shell profile 
around the Sun. 



Fig. 13 represents a comparison between the analytical 
solution and the numerical solution of the simplified BVP 
for a situation corresponding to a spherical sun and with 
a cosmological density of 10~^'^ kg /m^ for a model giving 
a thin-shell. We can see that the analytical solution is 
very good far from the sun but this agreement is less good 
inside the sun. Once again, for the determination of the 
PPN parameters, we only need the outside solution. The 
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analysis of Fig. 13 shows that the analytical solution can 
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Newtonian gauge [Xl where the coordinates are isotropic. 
More precisely, the PN parameters are defined with the 
Jordan-frame metric in isotropic coordinates as 




(64) 



where x is the Jordan frame isotropic radial coordinate. 
In all this section, we consider weak gravitational field 
and therefore we use a Post-Newtonian expansion. In 
this context, the symbol « emphasizes the fact that the 
equality is valid only to first Post-Newtonian order. In- 
stead of using isotropic coordinates, one can perform 
the transformation to express the space-time metric in 
Schwarzschild coordinates (this transformation for con- 
stants PN parameters is derived in [46] ) using the trans- 
formation 



FIG. 13. Representation of the evolution of the scalar field 
around the Sun for a model characterized by the parameters: 
A = 10"^ GeV, a = 3.3, fc = 1 and p^o = 10-^''fcg/m^ With 
these parameters, we obtain a thin-shell around the sun. On 
top: comparison of the field profile around the neighborhood 
of the Sun obtained numerically by solving the BVP problem 
related to the scalar field and by using the analytical solu- 
tion. On bottom: relative difl^erence between the numerical 
solution and the analytical approximation. One can see that 
the analytical solution is very good far from the sun but is 
less efficient in the Sun. 



D. Post-Newtonian parameters 

In order to constrain the parameters characterizing the 
theory with solar system experiments, it is useful to de- 
rive the post-newtonian parameters. As indicated in re- 
lation (52), the post-newtonian parameters for a tensor- 



scalar theory without potential is given by 

fc2 



7 - 1 = -2 



4Tr + fc2 ' 



^-1 = 0- 



As shown in the previous section, when a potential is 
added to the theory, the chameleon effect can reduced 
the coupling constant fc to an effective coupling constant 
fcoff as indicated in ( 63 ) . Considering the chameleon the- 



ory as a Brans-Dicke theory with an effective coupling 
constant fces as it was done in the original Chameleon 
papers [81 191 is not correct. Instead, one has to com- 
pute more carefully the PPN parameters following the 
procedure described in Damour and Esposito-Farese [5]. 
This procedure is used in Faulkner et al. [13] but they 
used a definition of the Post-Newtonian Parameters in 
Schwarzschild coordinates which can be misleading since 
the original PN parameters are defined in Standard Post 



~^~\^^ ^^^^ 



(65) 



A simple calculation gives 



I 1-2^ + 2(^-7) 
c r 




Pd^]^(66) 



where 7' ~ dj/df is the derivative of the PN parameter 
7 with respect to the Schwarzschild radial coordinates. 
This simple calculation shows that when we are dealing 
with space-time dependent PN parameters, one has to be 
very careful with the gauge in which the PN parameters 
are defined since its spatial derivative can appear. For 
example, in Faulkner et al. [13], they used Schwarzschild 
coordinates to identify a space dependent 7 but without 
any term related to 7'. This means that the expression 
they derive will differ from the one defined in isotropic 
coordinates which is the definition of 7 given in Will [T]. 

In the following, we will follow the procedure described 
in Damour and Esposito-Farese [S] being very careful to 
the gauge used. In the following, x will refer to isotropic 
radial coordinate and r to Schwarzschild radial coordi- 
nate and the ~ always refers to Jordan frame quantities. 
We start from the spherically symmetric metric in the 
Einstein frame in the weak field limit (in isotropic coor- 
dinates) : 



ds-" 



, Mb 
"mix 

, Mb 

'mix 



Mb ' 

mix, 



dt^ 



[dx' + x'dn'). 



(67) 



18 



Then, we can perform the conformal transformation to 
derive the Jordan frame metric 



ds 



^2 



(68) 



To transform the last expression into the form (64 1, 



one uses the field profile (61 1 as function of the 



Schwarzschild radial coordinates and the transformation 
between Schwarzschild and isotropic coordinates (r « x+ 
Mb/rrip) to express the field profile as function of x- Fi- 
nally, using the expansion A'^{(j}) — A^((/)oo)e^'^('^~'^=°^ = 
A^(l + 2k{(j} - 0oo) + 2fc2(0 - (/)oo)^) and the rela- 
tions between Jordan and Einstein frame quantities 
(i'' — A{(f)oa)x'^ — AoaX^ and mt, = Aoorhh), one obtains 
the metric (|64l with 



7 



cff 



47r 



„-moo(x-Xi.) 



47r- fcfcoffe-^-^x-X") 
Att + fcA:effe-™~(x-X(,) 

^ = 0- 



(69a) 

(69b) 
(69c) 



In practice, moo is tiny and the exponential decay can be 
neglected leading to the expression 



7 



Air — kkgff 
47r + fcfcoff 



(70) 



This approach gives the PN parameters as defined by 
Will [1, and avoid any confusion due to the choice of coor- 
dinates to define the 7 parameter. However, in a lot of pa- 
pers [21 [131 m] this derivation is done using Schwarzschild 
coordinates without the term linked to the derivative 7'. 
For this reason, we expect differences in the expression of 
the PN parameter 7 due to the difference of gauge used. 
This difference is present in relation (5.50) of De Felice 
and Tsujikawa [51 where an additional term is present 
in the denominator of 7. Nevertheless in [121 HI] ™oo is 
supposed to be very small and they neglect terms that 
are in fact closely related to the 7' contribution and they 
find a result similar to ( 70 1 (with some differences due to 
the choice of the starting action, see footnote [2 ) . Finally 
considering the theory equivalent to Brans-Dicke (which 
is done in Khoury and Weltman [SI [5] ) with a coupling 
constant fccff leads to a 7 = (47r — fc^) /{Air + fc^g )) which 
is quite different than expression ( |70| ) and which leads to 
a misevaluation of the PN parameters. 

The thin-shell effect has the property to hide the devia- 
tions from GR on solar system scales. Indeed, in the case 
of a thin-shell profile, the deviation from GR is quantified 
by the quantity 



7 - 1 = -2 



fcfcoff 



47r -|- fcfcoff 



-6 



47r + 3efc2 



(71) 



This shows that big values of the coupling constant k can 
be compatible with solar system experiments such as the 
measurement of the Shapiro delay by Cassini (531 if the 



k are completely excluded in Brans-Dicke theory. Using 
(62), the last expression becomes 



7 - 1 = -2 



fc('/>oo - (jjb) 
+ fc(0oo - (t>b) 



(72) 



It is possible to push the analytical development fur- 
ther by considering that (/)oo >> ipb (which is justified 
since it is shown in Tsujikawa et al. [15] that (poo /4>b = 
(Pb/poo)^^^"^^' - 1030/(«+i) if we use the cosmological 
and the sun density). Isolating / 
tion 



from the last rela- 



^b(7-l) 
2 + (7-l)^ 



(73) 



and using the Cassini constraint ( |53[ ) and the value of the 
Newton potential at sun radius ${, = 2.12 10~^, we get 
some boundary value fc^oo < 2.12 10~^^. The value of 
(poo is given by the value of the field minimizing the effec- 
tive potential VeS and is given by relation (37). For very 
low value of (j)oo (which is the case here), we can use the 
Taylor expansion of the W-Lambert function Wl {x) « x 
and we find the constraint 



aA4 



< 2.12 10" 



(74) 



Let us summarize the discussion about Post- 
Newtonian parameters. The 7 post-newtonian param- 
eters has been computed with great care of the gauge 
used. The 7 PN parameter is given by expression (70) 
where k^s depends on the thin-shell parameters e. If 
we have a thick-shell and the the- 



ory is equivalent to Brans-Dicke and the Cassini con- 
straint gives the usual constraint on the coupling con- 



stant k^ < 10" 



we have a thin-shell 



and the Cassini constraint is given by the constraint ( 74 ) . 

For the models within the 95% confidence region of the 
supernovae la likelihood analysis (see Section IIIC), we 
compute the thin-shell parameter e given by relation ( 62 1 . 
Fig. [14] shows the evolution of £ for model considered in 
Section [nrcj We can see that all these models present a 
huge value of the thin-shell parameter. This means that 
no thin-shell mechanism is playing any role for all these 
models and they are all equivalent to Brans-Dicke theory. 
Therefore, for models explaining cosmic acceleration, no 
chameleon effect is present on solar system scale. Since 
the theory is equivalent to Brans-Dicke, the constraint 
on 7 gives fc^ < 10~^. 

Another way to see the incompatibility between cosmo- 
logical observations and solar system observations is to 
take the (a. A) plane already presented on Fig. [3) With 
relations (l70| and (63), it is possible to compute the 7 



post-newtonian parameters for different models charac- 
terized by the three parameters A, a and k (considering 
a cosmological observed matter density of 10~^^ kg/m^) 



Fig. 15 represents the area in the parameter space (a. A) 



thin-shell parameter e is small enough while big values of that are compatible with 7 constraint ( 53 ) for different 
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FIG. 14. Representation of the thin-shell parameters e. The 
different lines represent different values of the coupling con- 
stant and the filled areas represent different values of a be- 
tween 0.5 and 3. 



values of the coupling constant. Several remarks need to 
be done. First of all, we see that high coupling constants 
are not completely rejected by solar system experiments. 
For example, fc = 1 or fc = 10"* are completely rejected 
by considering massless tensor-scalar theory while it can 
be seen that these coupling constants can be admissible 
with respect to post-newtonian test of gravity in a certain 
domain of the space parameters (a. A). In particular, for 
low values of A and high values of a, the chameleon ef- 
fect is strong enough to reduce the scalar charge kes so 
that the 7 parameter (70) satisfies the Cassini constraint. 
This was the breakthrough of the original papers [HI [H] ■ 



Nevertheless, we can see on Fig. [15] that it is not possi- 
ble to explain cosmic acceleration and to pass the solar 
system experiments at the same time (at the exception 
of the small coupling constants k which is also the case 
for Brans-Dicke theory). Indeed, there is no intersection 
between the curves representing the relation between A 
and a imposed by the cosmological observations and the 
domain allowed by the PPN constraints. 

In the calculation of the PPN parameters in Fig. [15] 
the density at infinity of the solar system is assumed to 
be the cosmological density {p{r = 00) = po). This situ- 
ation is a simplified situation. In reality, the situation is 
much more complex since the solar system is not plunged 
into a cosmological fluid but is a part of a galaxy which 
itself is a part of a galaxy cluster and so one. . . As a 
consequence, the scalar field will have an evolution at 
all the different scales from solar system scales to cos- 
mological scales (passing by the galactic scale and the 
galaxy cluster scale and so one. . . ). All this modeling is 
quite complex and beyond the scope of this paper. Nev- 
ertheless to have some intuition, we compute the PPN 
parameter by assuming that the density at infinity of the 
solar system is given by the galactic density at a dis- 
tance of the solar system. We assume a galactic density 
of Pgai ~ 10~^^kg/m^ at solar system distance from the 
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FIG. 15. Representat ion of couples (a, A) admissible with 
the Cassini constraint ( 53 1 for different values of the coupling 
constant k and representation of the couples (a, A) needed to 
explain cosmological observations. The density at infinity of 
the solar system p{r = 00) = 10~^^ kg/m'^ which corresponds 



to the cosmological matter density. 



galactic center. Fig. [T6| represents the admissible domain 
in the (a, A) plane supposing that the density at infinity 
is given by the galactic density p(r = 00) = Pgai- We 
can see that the situation is better than considering the 
cosmological density. The area of the region satisfying 
the PPN constraint are bigger but there is still no inter- 
section between the cosmological relation and the solar 
system constraint. 
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FIG. 16. Representation of couples (a. A) admissible with 
the Cassini constraint ( 53 1 for different values of the coupling 
constant k and representation of the couples (a, A) needed 
to explain cosmological observations. The density at infinity 
of the solar system p{r = 00) = \Q~^^kg/rn'^ which roughly 
corresponds to the galactic matter density. 



The results of this analysis is that, on solar system 
scales, the chameleon effect can allow high coupling con- 
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stant k that are normally excluded by Brans-Dicke anal- 
ysis at the condition that parameters a and A are in 
the admissible domain (which mean for high a and small 
A). But we have shown that there are no intersection 
between models explaining cosmological expansion and 
models authorized by solar system experiments. 



V. CONCLUSION 

In this paper we have studied massive tensor-scalar 
theories at two different scales: on cosmological scales 
and on solar system scales. We have presented combined 
constraints of the so-called chameleon field introduced by 
Khoury and Weltman |H1 IS] ■ In this communication we 
have only focused on the original model which is charac- 
terized by an exponential coupling function A{(j)) = e'^'^ 
and on a Ratra-Peebles [20] run-away potential. We have 
also derived a unambiguous definition of observables like 
luminous distance and 7 Post-Newtonian parameter, 
leading to different predictions than in literature. 

On cosmological scales, we have derived the cos- 
mological equations of evolution by using the Jordan 
frame density parameter as source term in the Einstein 
equations. Using the latest Supernovae la data [30], we 
have derived confidence regions for the coupling constant 
k, for the parameters characterizing the potential and for 
the cosmological matter density today ilmo- Considering 
models within this confidence region, we have derived 
a relationship between the two parameters character- 
izing the Ratra-Peebles potential. The relationship 
found ( 33 1 is exactly the same as the one obtained 



in quintessence model (characterized by a vanishing 
coupling constant fc = 0) [32]. A detailed analysis 
of the cosmological evolution has been presented. In 
particular, the contribution of the scalar field has been 
parametrized by an effective fluid of dark energy in GR. 
We have presented the evolution of the effective equation 
of state parameter lode and we have shown that the 
derivative of this parameter can become positive which 
can lead to a test to distinguish tensor-scalar theories 
from quintessence scenarios. Moreover, we have shown 
that for high coupling constants, the cosmic acceleration 
is mainly produced by the non-minimal coupling and 
not by the scalar self interaction. Finally, we identified 
two dynamical regimes that can explain the cosmic 
acceleration. This two regimes are characterized by two 
minima in the curves as shown in Fig. [2] 

In the second part of the article, we have studied the 
solar system configuration modeled by a spherical and 
static central body (the Sun). We have computed nu- 
merically the field profile solving the full boundary value 
problem composed by the Einstein field equations and we 
have compared these numerical simulations to analytical 
solutions that can be found in the chameleon literature. 
In particular, we have shown that the analytical solu- 



tions are very good outside the central body but are not 
completely satisfactory inside the body. These simula- 
tions have confirmed that the thin-shell mechanism can 
appear in the surrounding of the Sun. Using analytical 
expression outside the central body, we have derived the 
post-newtonian parameters. In this derivation we have 
taken a great care to the gauge used to define PN param- 
eters comparing our approach with existing expressions 
in chameleon literature. In the case of a thin-shell sit- 
uation we have confirmed that the scalar charge of the 
central body can drastically be reduced by the relation 
fcoff = 3eA: where e is the thin-shell parameter. Contrary 
to what is claimed in Khoury and Weltman [HUH], we have 
shown that the theory is not equivalent to Brans-Dicke 
theory with an effective coupling constant fcog. As a con- 
sequence, the 7 constraint is different from what was de- 
rived in the original chameleons papers. Nevertheless, we 
have confirmed that in the case of a thin-shell situation, 
the reduction of the scalar charge implies a reduction of 
the deviation of the post-newtonian parameter 7 from 
its GR value. This confirms the fact that high values 
of the coupling constant k can eventually pass solar sys- 
tem constraints while these coupling constants were pre- 
viously excluded by considering traditional tensor-scalar 
theories. 

For models within the 2(t cosmological confidence 
region we have computed the thin-shell parameters e 
and we have found that this parameter is always very 
huge compared to the unity (e >> 1). This means 
that no thin-shell effect can be invoked on solar system 
scales for models fitting Supernovae la data and as a 
consequence no reduction of the scalar charge can be 
observed. The theory is therefore equivalent to Brans- 
Dicke and the Post-Newtonian constraint implies that 
only small coupling constants are admissible. Another 
presentation of the combined constraints have been 
presented in the form of admissible region in the (a. A) 
plane (where A and a are the parameters characterizing 
the Ratra-Peebles potential). In this plane, we have 
shown that there is no intersection between cosmological 
admissible domain and solar system admissible domain. 

The main conclusion of this work is that the chameleon 
field (if one consider the initial model by Khoury and 
Weltman [H |9] ) can not explain cosmological accelerated 
expansion and pass solar system constraints at the same 
time. On one hand, we have shown that there are models 
fitting Supernovae la data relying most on non minimal 
coupling than on self interaction potential. On the other 
hand, there are models satisfying solar system constraints 
even for high coupling constant (which are rejected if no 
chameleon mechanism is present). But there are no mod- 
els satisfying both conditions. We emphasize that this 
conclusion is only valid for exponential coupling constant 
and for a Ratra-Peebles potential. It would be interest- 
ing in a future work to study other type of potential (like 
the SUGRA model [47]) and other type of coupling func- 



tion (like the gaussian coupling constant A((j)) 
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which can eventually lead to a different conclusion. 
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